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ABSTRACT 


This paper proposes and evaluates a general search 
functicn which can be used for determining the 
erficiency of a search in which the searcher is 
required to visit the maximum number of targets in an 
qimored period of time. A lower bound on this 
Tuncticn is presented as a consequence One the 


“Semimyopic" strategy utilized in the Simulation. 
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I. INTRODUCTION 


Two of the factors that must be considered in planning 
the search of a large area are the search radius (maximun 
detection range) and search speed. Due to a decrease in the 
Maximum detection range of some sensors with increasing 
search speed, the interaction of these two factors and their 
ultimate erfect on search efficiency must be accounted for. 
The development of some analyticai means by which the 
Planmuer could account for the effect of changing search 


radius on search efficiency prompted this research. 


The concept investigated was that of a single searcher 
reguired to visit aS many targets as possible, in an 
infinite field of uniformly distributed targets (poisson 
field), in an alloted period of time. The searcher was 
placed randcmly in the field of targets at the start of the 
search. Visiting the same target more than once was not 
allowed. Search speed was assumed to be much greater than 
target speed and, hence, all targets were considered 


stationary in the field. 


The search strategy employed becomes extremely impcrtant 
as the séarch radius increases with respect to expected > 
target spacing. It seems logical to always go to the nearest 
detected, but unvisited, target; however, such a strategy 
allows the searcher to backtrack and travel over previously 
covered areas. It may weli be better for the searcher to go 


tOoewa Lresh area. 


Figures 2 through 7 snow the typical track of a searcher 


employing a "go to the nearest target" or “semimyopic" 





Strategy for increasing values of search radius. In each 
case the searcher started from the center of a poisson field 
of 1000 targets. If no targets were within the search 
radius then a search course was taken away from the center 
of mass of all previously visited targets, until a target 


was detected. 





If. BACKGROUND 


The objective of the searcher is to maximize the number 
of targets visited in a unit of time, subject tc the 


following assumptions: 
ee Soecaren f1ceid is of infinite size. 


spibe = target] are Uniformly distributed over the 


e The targets are stationary. 


This search froblem can be divided into essentially three 


G@ieetanet Situations or cases. 


Peer = SEARCH RADIUS MUCH LESS THAN EXPECTED TARGET 
SPACING 


In the limit this is the classical random search problen 
@escriped in Ref. 1. If the search radius (r) is very 
Small, then the number of targets visited is approximately 


equal to the number of targets detected. 


Pewpcaoe £2 — StaAkRCH RADIUS NUCH GREATER THAN EXPECTED 
PAMGET SPACING 


In the limit this approaches a situation of totally 
unrestricted visibility, and the searcher can pick and 


Cheese nis optimum track. 





a 





See coon. Ili = 


This is 


ANY OTHER SEARCH RADIUS 


by 


far the most interesting case and the 


cbhjective implies finding an optimal search strategy when 


anything more 


sophisticated then random search can be 
employed. Possible search strategies could include any one 
om, Of 2 combination of, the following: 


e Random searcn. 


» Search along a steady course. 


a Patterned search. 


» GO 


es GC 


to 


to 


the nearest target if one or more targets 


the nearest area of "high target density" 


1£ more tken one held. 
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Pi NA LY SES AND THEORY 


a ae ee ee ee ee ee ewe oD eee ee eee oe 


If the number of targets visited is fixed, then the 
average tctal distance traveled {1) to visit all of the 
targets is a function of the number of targets visited (n), 


the search radius (r), and the target density (p). Thus, 


(Ge oie ie Aue 


Ae DIMENSICNAL ANALYSIS 


By digensional analysis, 
BL Ga) ee icw/S) lat) 


since all distances can be measured in terms of 1/Vvp. 


Assuming that 1(r,p,n) 1S asympotically linear in n, then 


lim’ 2£ (cVp,n) /n = 1/g (rvp) 


n>@w 


The primary cbjective of this thesis was the discovery of 


this dimensicnless function, g(x), where x = Ivp. 


Paee kPECTATICNS FOR g(X) 


ec y 


Se th 
C= Distance From (1-1) TOs Target 
1 


Mnere 1 =i 2, 3) ss 
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Then, 


P(D 2d) = P(No fTgts in Trdée) 
Siar p 
= 1-e 
and, 
co 
BG 
Ete ={e dy = 1/(2vp) p> 0 
° 


Since the first detection is based on the searcher being 
Piesecdsan a VirGin field, it follows that the probability of 
being clcse to an unvisited target is greater at the start 


then at any other time of the search. Thus, intuitively, 


Eee Se (De) lee ee 
1 a 
Buty 
i= 2. Decree De) 
1 n 
=O 
IWigt c=) ene oD eet De) y p/n 
n?co 1 n 
n 
= lim 2.35 2(D.)Vp/n 
neo i= aL 
> lim 2nE(D )Vp/n = 2E(D )Vp = 1 
nh »c 1 1 
Thus, 


Gx os 


meuexye is small, then ome can conclude that classical 
random search, as described by Koopman [Ref. 1], applies and 
that the number of targets visited is approximately equal to 


the number of targets detected. Thus, applying the standard 


2 





“cookie cutter" model, gives 
N = ap = 2rip = 21Vp(rvp) 
Where a= Area Searched 


and for sail x, 


Sawa BSOTIMATOR FOR g(x) 


bet, 
j oe i OM ee 
then, 
2) (A es et t 22020 + Been 
Soe ay ry a ae 
a 1/g (Xx) tien Ls Large 


Thus, for large n, Z is an eStimator for 1/g(x). 


This thesis used 1/2 as an eStimator for g(x) even 
though E(1/Z) 2 g(x)’. Because the Var(Z) was small when ao 


was large, 


BAA7 2 yee EZ) =. G(X) 
The actual estimator utilized was, 
g(x) = N/(2LVp) 


Where N = Number Of Targets Visited 


en 
Total Distance To n Target 


L 
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The development of optimal search strategies for the 
full spectrum of search radii to be employed was dismissed 
frcem the outset, and two practical search strategies settled 
On. The resulting strategy may best be defined as being a 
semimyopic strategy. This semimyopic strategy provided a 
lower bcund on g(x), and hence on the expected number of 


targets that can be visited by a searcher. 


pee CO DERECT MOWE STRAT ECY 


— 


If there are targets within the search radius of the 


searcher, he visits the nearest target. 


eee ero TReCT SYOVE STRATEGY 


ae no targets are held initially, the search is 
conducted along a randomly determined course. If, having 
visited at least two targets, there are no targets within 
the search radius of the searcher, the search course is 
determined ry the position of the last target visited and 
the center or mass of all previously visited targets. In 
both cases, the search course is always directed towards 
areas not previously searched and the searcher always visits 


the first target detected. 


14 





Pigure 1 shows the plot of g(x) vs x and its one-Sigma 


boundaries oktained by computer simulation. 


Assuming that a search is to be conducted in a poisson 
Mier Of targets, the planner can, by uSing the relaticnship 
developed earlier and Figure 1, determine the expected 
number of targets that can be visited by passing close 
akcard. Furthermore, by knowing the effects that search 
speed nas on the search radius of the sensor to be employed, 
he can ccmpare expected search results and select the most 


efficient or favorable speed at which to conduct the search. 
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g (x) 


#igUre Ine Sgitx) VS X 
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VI. SIMULATION MODEL 


A unit square with 10,000 uniformly distributed randon 
targets was used to represent an infinite field, and a 
Simulation séarch program was written in FORTRAN and run on 
the IBM 360. The searcher commenced his search from the 
center of the unit square and, after a series of direct 
(target within search radius of searcher) and indirect (no 
target within search radius of searcher) moves, would 
terminate when (1) the searcher was within one search radius 
or the fields boundary, or (2) the distance to the next 
target was greater then the distance to the boundary of the 
field. Appendix A contains a flow chart of the models basic 


POmEEC . 


A minimum of 40 runs were made through the same randon 
fields for various values of search radius, eguivalent to x 
ranging Beene O.025 tO <Lntingey. The results of the 
Simulaticn runs are contained in Appendix B, and Figure 1 is 


tPicmi@esiulting plot of g(x) vs x. 


The effect of x, and consegquently search radius, on the 
searchers track and the number of targets visited are shown 
mietiauee= 2 through 7. The 1000 target field 15 identical 


Mma each case. "Direct" refers to the Direct Move Strategy. 


The results from all runs made in 10,000 target fields 
miercated that for x 2 1.0 the Direct Move Strategy 
prevailed, and that for x 2 4.0 the searcher always held the 


next target. 


a =F 





The simulation required 170K on the IBM 360 and the tine 
reguired for a run through 40 poisson fields of 10,000 
tarcets each varied from about one minute, for small values 
of search radius, up to about fifteen minutes, for large 


values of search radius. 


13 





VII. CONCLUSION 


aS oe oe ee ee 


The results obtained for small values of x correlate 
nicely with random search theory. Based on random search 
theory, it was expected that g(x) would egual x for small 
values of x, and Figure 1 shows that this relationship holds 


@orevalues Of X Up tO approximately 0.5. 


It was anticipated that g(x) would approach some fixed 
value with increasing x, yet Figure 1 and the data obtained 
indicate a slight peak at about x = 2.0. The only plausible 
explanation for this peak is that the searcher is penalized 
for emplcying a myopic strategy in a hypermetropic 


Situation. 


A close inspection of Figures 2 through 7 leads one to 
conclude that a strategy of going to the nearest area of 
high target density could substantially increase the number 
of targets visited per unit distance traveled, and 
Subseguently g(x). To verify this, seven visual trials were 
made by three individuals using the same 1000 target field 
used in Pigures 2 through 7. The objective of the 
individuals was to visit as many targets aS possibile with 
unrestricted visibility. The resulting values of g(x) ranged 
Pao. 7622 to O.8401 with g(x) = 0.8075. By computer 
Simulation, for unrestricted visibility, g(x) = 0.6720652 
Tor this field. These results tend to substantiate the 
belief that a searcher employing a myopic strategy in a 
hypermetrcpic situation 1s penaiized, and the earlier clain 
that the semimyopic strategy served to provide a iower bound 


On g(X). 
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APPENDIX A 


FLOW CHART FOR SIMULATION MODEL 


SET: SEARCH RADIUS (RW) 
RANDOM NR SEED 
TGT DENSITY (N) 





INITIALIZE 


GENERATE FIELD 
OF NHN iGis 


DETERMINE TG6Ts 
IN SUB VAR AL OE 
RADIUS RA 
FROM SEARCHER 


ES DAE Roe a Gl 
Wi tae 8 OF 
S EARCHER? 








(DIRECT MOVE) 


AREA EOUNDARY 


mS SUE 
Groner LHAN LGL? 


NO 







BOUNDARY 


IS NEAREST BOUNDARY 
CLOSER THAN wWEXT 
TARGET? 






YES 


155 
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UPDATE SEAR 
POSIT IN F 


t 


Deora IN Sep es TANCE 
FROM CENTER OF SUE AREA 
TO NEAREST BOUNDARY 


CHERS 
cae} 







GO TO TGT 


UPDATE SEARCHERS 
POSIT IN SUB AREA 


RECORD DATA: 
PCR EAENY N 


DISTANCE TRAV 








(INDIRECT MOVE) 


TS oo FARCH Oars 
RW OF 


ER 
UB 
DA 


DETERMINE DISTANCE 
FROMSGENT ER SOL eSUs 
AREA TO NEAREST BOUNDARY 





NO 
SE DETERMINE DISTANCE 
FROM SEARCHER TO 
NEAREST BOUNDARY 
NO IS NEAREST BOUNDARY 
3 a THe Ne RWeOb es EARCHUER? 


is) 


D@bSmoUB AREA CONTALEN | NO 
NEAREST BOUNDARY? GO TO TGT 
UPDATE SH#ARCHERS 
POSIT PNesUs AR EA 


RECORD DATA: 
INCREMENT NR VISITED 
DISTANCE TRAVELED 







DEER ENE DiS TANCE 
FACM NEXT GT 
TO NEAREST BOUNDARY 










Poel obane s 50 
NEAREST BOUNDARY 
ios LHAN: Re? 






EES 


SD Neg 









COMPUTE g(x) 


COMPUTE (sta bea tucs 
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APPENDIX 5 


COMPUTER OUTPUT 


Number g (X) Sigma 
Of Runs See eae yea 
40 OnOm2uS 5. 020315156 
40 0.0678127 0.0360265 
40 0.0753077 0.0260168 
40 0.1006204 0.0282282 
80 0.1318160 0.0342952 
80 0.2574581 0.0476352 
80 0.4808525 0.0567625 
80 0.6374704 0.0501638 
80 0.6741247 0.0418967 
80 0.6959368 0.0319426 
80 0.7058284 0.0312344 
80 0.7143928 0.0279759 
80 0.7196197 0.0244312 
80 0.7194164 0.0238991 
80 0.7176019 0.0225489 
40 0.7181079 0.0200270 
40 0.7136417 '0.0247005 
40 0.7161400 0.0227902 
40 0.7124111 0.0252671 
40 0.7110950 0.0260190 
4.0 0.7114341 0.0255649 
B0 0.7117410 0.023186 1 


TABLE I 
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